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We present our results on meson and nucleon screening masses in finite temperature two flavour 
QCD using smeared staggered valence quarks and staggered thin-link sea quarks with different 
lattice spacings and quark masses. We investigate optimization of smearing by observing its effects 
on the infrared (IR) and ultraviolet (UV) components of gluon and quark fields. The application of 
smearing to screening at finite temperature also provides a transparent window into the mechanism 
of the interplay of smearing and chiral symmetry. The improved hadronic operators show that above 
- - ■ the finite temperature cross over, Tc, screening masses are consistent with weak-coupling predictions. 

' There is also evidence for a rapid opening up of a spectral gap of the Dirac operator immediately 

! above Tc. 

O . 
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I. INTRODUCTION 



Screening masses control finite volume effects at finite temperature in equilibrium. Studies of the final state of 
d ' fireballs produced in heavy-ion collisions indicate that they are near equiUbrium. So the study of screening masses a 
little below the QCD cross over temperature, near the freeze out, should improve our understanding of experimental 
Qh' conditions. In addition, the vector screening masses below Tc should be of direct relevance to the study of mass 
^ [ spectra of dileptons and photons. 

There are also interesting questions about the nature of the high temperature phase which are addressed by a study 
of screening masses. In QCD at temperatures, T, of a few times the crossover temperature, analysis of the weak 
coupling series in powers of the gauge coupling, indicates that the physics of the magnetic scale of momentum, g^T, 
K*" ■ is potentially non-perturbative. As a result, it may be possible to find phenomena in hot QCD, only involving harder 
scales, which are amenable to a suitable weak coupling analysis. For example, the fermionic part of the pressure, as 
well as its derivatives with respect to chemical potentials, the quark number susceptibilities (QNS), seem to admit 
reasonably accurate weak-coupling descriptions at temperatures of 2Tc or above [l| . 

However, even among static fermionic quantities, screening masses (the inverses of screening lengths) present a 
^ii^ ■ confused picture. Most computations have been performed with staggered quarks, and these seem to indicate that 
, there are strong deviations from weak coupling prediction ■ On the other hand, computations with Wilson quarks 
give results which are closer to free field theory , although they deviate in detail from predictions of weak coupling 
theory [1,0. Since the same pattern is visible in the quenched theory [lo| . we can attribute the major part of the 
discrepancy to valence quark artifacts. Here we examine this question systematically using staggered sea quarks and 
improved staggered valence quarks. Indeed, we see that smeared valence quarks provide a significant improvement. 
• Using these we find that a weak coupling expansion does work quantitatively for the description of fermionic screening 
masses at finite temperature. In addition, our results may constrain models of thermal effects on hadrons below and 
close to the QCD cross over. 

A significant technical component of this work is the exploration of the cause of im pro vement in lattice measurements 
when smeared gauge fields are introduced into the staggered quark propagators [Ill4l4| . Smeared operators have been 
explored extensively in the literature earlier [isj . Here we explore optimization of smearing parameters by direct 
observation of the effects on UV and IR modes separately. It also turns out that the application to finite temperature 
provides a transparent window into the interplay of improvement and chiral symmetry. 

Discussion of technical lattice issues in this paper are confined to the next two sections. Readers who are interested 
only in the results for thermal physics can read the last two sections. 
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TABLE I: The number of independent configurations, A'^, obtained with the couphng, /3, the bare quark mass, am, and the 
autocorrelation time, r, for that simulation. Also given are the plaquette value, P, measured at T = 0, and the temperature, 
T /Tc, inferred from it. 



II. METHODS AND DEFINITIONS 



We generated configurations for the Wilson gauge action and two flavours of thin-link staggered sea quarks using the 
R-algorithm. For am = 0.015 we used lattice sizes Nt x with TVt = 4 and Ng = 8, 12, 16 and 24 for finite T studies, 
and scanned a range of gauge couplings, /3, to find the cross-over coupling /3c. This is completely standard, and the 
results are collected in Appendix 1X1 The simulations were done using a MD time step dt = 0.01 and trajectories with 
number of steps, Nmd — 100(Afs/8). We checked that halving the time step did not change the results. We observed 
that it was sufficient to discard the first three hundred trajectories for thermalization. The configurations analyzed 
were thermally equilibrated and spaced one autocorrelation time apart. Details of the runs and statistics are collected 
in Table m This data set is called the set N in the rest of this paper. 

We also studied configurations generated earlier along a line of constant with Nt = A defined by setting 
am = 0.025 at the corresponding /?c = 5.2875 [iBl and Nt = 6 defined by the choice am = 0.01667 at its f3c = 5.425 
[l7| . The data from is referred to as set O in this paper, and the data of as set P. Hadronic screening masses 
from the data set P have been reported earlier using thin-link staggered valence quarks Q ; its inclusion in this study 
enables a clear understanding of the effects of smearing. 

We studied screening correlators of mesons and the nucleon. The valence quarks were improved using one level 
smeared gauge links |lll4l4j : the optimization of the smearing algorithm is discussed in Section Hill In the course of 
this study we needed to estimate the extremal eigenvalues of the staggered Dirac operator. This was done using a 
Lanczos iteration [lH . The tridiagonal matrix generated using this process was diagonalized using the Lapack routine 
DSTEVX. The investigation of smearing also needed the determination of the taste partners of the pion. For all 
the correlation functions we used Coulomb gauge fixed wall-sources to project on the modes with vanishing spatial 
momentum. At T = and for temperatures below Tc, multiple wall-sources separated by 4 lattice units were used. 
We checked that these gave statistically independent results; an observation that could be justified after the fact by 
the measurement of the Goldstone (local) pion mass. 

The screening correlator for the meson 7 was parametrized as 

Cj{z) = Ay cosh 

The alternating component is absent for the Goldstone pion [l9|. Among local operators, we measured the scalar (S) 
corresponding to the a/ao meson at T = 0, the pseudoscalar (PS) corresponding to the tt at T = 0, the vector (V, p 
at r = 0) and the axial vector (AV). At T = all three polarizations of the V and AV are equivalent. However, for 
T > 0; we need to distinguish between the spatial (Vg, AVg) and temporal (14, AVt) polarizations. For the study of 
taste symmetry we also measured the non-local taste partners in some of these channels. Following p^ . the nucleon 
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FIG. 1: Power suppression, Q, in the UV and IR with different kinds of smearing at three different lattice spacings for T — 0. 
A halving of the lattice spacing leads to a weak change in the optimal value of e. 



correlator is parametrized as 
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The screening masses, fj,j, fi^, and the remaining parameters were extracted from the measured correlators by fitting 
to the above forms. The covariance between the measurements at different z were taken care of in the fits. The 
mean and the error of the parameters were estimated by bootstrap. Fits were made to the ranges Zmin ^ z < Zmax, 
where Zmux was never more than two sites from the middle of the lattice, ^min was never less than two sites from the 
source, and the number of data points used was always greater than the number of parameters being fitted. Among 
the fits satisfying x^/DOF < 2, we chose as the reported estimate of the parameter and its error to be that which 
was consistent with the smallest fi within 2a and had the smallest error. Chiral symmetry restoration can be tested 
through the mass splittings 



A, 



fJ-s - fJ-PS, 



(3) 



as well as the parity projected correlators C^'' = ± (—1)^(7'*' , where the mesons 7 and 7' are parity partners 

We computed quantities in a fermionic free field theory (FFT) by numerical inversion of the fermion matrix on a 
trivial gauge configuration (all links being the unit matrix). These quark propagators were then subjected to exactly 
the same analysis as in the interacting theory. The negative chiral projections of the screening correlators, C~'^, 
vanish in the chirally symmetric phase, and the approach to FFT can be studied using the positive chiral projections. 



III. STUDY OF SMEARING 



Smeared gauge links improve the scaling behaviour of staggered quarks (15| . This has been measured through 
staggered pion taste splitting. Optimal parameter values have been obtained numerically and there have been attempts 
to understand the results in weak-coupling theory [l^ . 
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FIG. 2: Speedup of the conjugate gradient inversion for T = at two different lattice spacings as a function of e in the various 
different smearing schemes. The optimal value of e in each smearing scheme agrees with that seen in the glue sector. This 
closely follows the change in the condition number, k, of the fermion matrix. The last panel shows that the change in the 
condition number comes from the UV, i.e., Xmax', the IR, i.e., Xmin, is almost unchanged by smearing. 



We examined four schemes which are currently popular: APE [Tl|, HYP [T^, Stout [T^, and HEX All these 
schemes involve replacing the gauge field on a link by a weighted sum of gauge transporters over different paths 
connecting the end points of this link. The more steps of such smearing we take, the more non-local the action 
becomes. In order to retain a degree of locality compatible with the sea quark action, we restricted ourselves to 
one step of smearing. The APE and Stout schemes have a single free parameter, e, which determines how much 
importance is given to link neighbours. The HYP and HEX schemes have three different fattening parameters in 
three orthogonal directions. We restricted our study to the subset which have equal contributions from all directions, 
controlled by a single parameter e. 



A. Optimization of smearing parameters 



The usual lore about smearing is that it suppresses the dependence of operators on high-momentum field modes. 
Since the lattice cutoff affects UV modes strongly, the result could be closer to the continuum limit. Since field 
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TABLE H: The best e for two different a, the second being half of the first, evaluated in different schemes and by different 
optimization criteria. The optimum parameter value in each scheme is nearly independent of a. 



operators have a gauge dependence, it is hard to test this idea directly on gauge fields. Instead we tested it on the 
plaquette at a site x averaged over all 6 orientations, P{x). As for any local operator, one can work with the Fourier 
transform, P{k), and the power spectrum, E{k), where 

P(fc) = ^ exp (ifc • x) P(a;) and E{k) \P{k)\^ , (4) 

X 

the mode numbers k^ = 7r(2£p + (p,)/N^, is the size of the lattice in the direction /x, the integers < < A^^, 
and Cm = periodic boundary conditions and 1 for anti-periodic. Periodic or anti-periodic boundary conditions 
imply that the independent modes are those with inside the Brillouin hypercube whose body diagonal, BD, joins 
the corners (0,0,0,0) and Ny/2, Nt/2). 

We used this power spectrum to find how smearing affects the UV and IR modes. We separated the IR and UV 
using hyperplanes perpendicular to BD. All modes within the Brillouin zone closer to the origin than a hyperplane 
am were called IR modes; conversely all modes within the Brillouin zone closer to the far corner than the plane (t^v 
were called UV modes. Everything else was a generic mode- neither IR, nor UV. We defined the suppression of power 
in the IR and UV as a function of e 

where Euv{^) is the power summed over all modes in the UV for a fixed value of e, and -©/^(e) is a similar quantity 
obtained by summing over all modes in the IR. The definitions of IR and UV are arbitrary, so one needs to check 
whether the results are sensitive to this definition. We placed the planes ctjh and auv at a fraction d of the length of 
the diagonal (with < d < 0.5, so that no mode is simultaneously in the IR and UV) from the nearest corner, and 
varied d. We observed that results were insensitive to d. 

We investigated Q numerically with thermalized configurations at T = using (3 = 5.2746 and am ~ 0.015. Periodic 
boundary conditions were used so that all C^ = 0. The variation of Qx with e is shown in Figure [1] One sees that the 
slope of the curve for Quv always starts off larger than that for Qm. Also, the slope of the latter seems to be close 
to zero. This shows that smearing can be used to modify the UV without modifying the IR. One can use this to seek 
an optimum value of e, such that Quv is as small as possible. In a simulation with dynamical smeared quarks one 
would have to do this without making a significant change in Qir. In this study the smearing is quenched; the set of 
gauge configurations is not changed by smearing, only valence fermions are affected by smearing. So in this context 
we are free to drop the condition on Qi^. 

We also investigated the quark mass and lattice spacing dependence of Qir and Quv by studying thermalized 
configurations at T = using /3 = 5.2875 and am = 0.025 as well as (3 = 5.53 and am = 0.0125. The first set has 
almost the same lattice spacing as the one with (3 = 5.2746, but has a somewhat different pion mass. The last two 
sets have the same pion mass but lattice spacings which differ by a factor of two. We show the results in Figure [T] 
As can be seen very clearly, there is a change in the overall suppression of power in the IR and UV, but the change 
in the optimum e is not large even when the lattice spacing is halved. The optimum values of e move down slightly. 
This movement is compatible with the intuition that finer lattices require less improvement. 

Interestingly, conjugate gradient inversion is also optimized at similar values of e [20| . In Figure [5] we show the 
number of conjugate gradient iterations required to invert a smeared staggered Dirac operator, Ncg, in a representative 
configuration drawn from thermalized ensembles. These results were obtained with a CG stopping criterion that the 
norm of the residual is less than 10~^^/V. Note that the performance of the APE and HYP smeared operators are 
very similar to each other, just as before. The behaviour of the Stout and HEX smearing are also similar, but quite 
different from the previous pair. Again, the lattice spacing and pion mass seems to make little difference to the 
optimization. 



6 



Using the smeared staggered Dirac operator, D, we found the minimum and maximum eigenvalues of D: \min 
and Xmax- We defined the condition number K(e) = Xmaxi.^) / ^min{^) ■ One expects that the number of CG iterations 
is closely related to K(e), as indeed it is seen to be (see Figure [2]). We found that Ami„(e) is independent of e to 
better than 1%, as expected, so Aniin(e)/AmjTi(0) is flat. The dependence of k on e is essentially due to the variation 

of Amax(e)- 

We can use any of these criteria, namely, the minimization of QuVi Ncg, or Xmax, to choose the best value of e. The 
results are shown in Table |lT] in different smearing schemes at different lattice spacings. There is a marginal decrease 
in the best smearing parameter in each scheme with decrease in lattice spacing. Wc sec that there is reasonable 
agreement between the best values obtained through the three methods. Given this, we choose to work with the 
values e = 0.6 for APE and HYP, and with e = 0.15 for HEX and e = 0.1 for Stout. 



B. Smeared quarks and chirally symmetric correlators 
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FIG. 3: Screening correlators from data set N above and below Tc- The signs of chiral symmetry restoration are clear with 
either thin- link or improved valence quarks in the form of pairwise degeneracies of correlators above Tc- However, improved 
correlators show even higher degeneracies at high temperature. Similar results are obtained for data sets O and P. 

Chiral symmetry restoration in the high temperature phase of QCD is easily seen in hadronic correlation functions. 
Below Tc the local meson correlators: S, PS, V, and AV arc quite distinct, but above Tc they collapse into one 
(see Figure [3]). A pairwise degeneracy of the S/PS and V/AV shows chiral symmetry restoration, and has been 
demonstrated earlier as well with thin-link staggered valence quarks. However, the near-degeneracy of the two pairs 
at high temperature, visible only after smearing, is a new observation. This occurs in all the data sets: N, O, and P. 

Pairwise degeneracy arising from chiral symmetry restoration is most easily seen in the vanishing of C^"'^-', C*^"^'^ 
and C*-~^=^ at high temperature On examining these combinations, it turns out that the degeneracy for T > Tc 
becomes clearer with smearing. For example, C''~^*'(z — 1/T) is (6 ± 8) x 10~^ at Tc with thin-link valence quarks, 
but becomes (0±2) x 10"'^ when optimal HYP smeared valence quarks are used. The improvement is most remarkable 
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FIG. 4: The correlator C' ^""^ shows interesting short distance (2 < l/T) spatial structure slightly above Tc (left). This effect 
barely persists into the hot phase (right). The results are displayed for set N. Sets O and P show similar behaviour. 
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FIG. 5: The ratio of chiral projections 0^+^"'' (left) and C'"'"'^' (right) at 2Tc to the respective FFT prediction with data sets 
O and P. The smeared correlators come close to FFT in both cases, whereas the unsmeared C^'^^^ is quite different. 



in the S/PS sector, where we found C^'^^z = l/T) = -3.3±0.1 at using thin-link valence, but -0.57±0.04 using 
optimal HYP smeared valence. At larger T all the negative chiral projections vanished. It was seen earlier Q that 
fQj. X > Tc vanished when z > l/T, but remained non-zero at short distances. In Figure|4]we show this effect 
at Tc and also, that it vanishes at 2Tc. A more detailed view of the temperature dependence is exhibited by showing 
how C*^~^=^(z = l/T) changes with T. Below Tc the correlator does not vanish, but the spatial structure seems to 
have entirely disappeared for T > l.OSTc. This gives one definition of the width of the chiral crossover; it is larger 
than the one implied by A/3c (see Appendix [A]) . 

(j(+ys) jg close to the FFT prediction with either thin-link or smeared valence quarks. With thin-link staggered 
valence quarks, we see that C^~^^^ is different from FFT, as previously observed. However, on smearing, they become 
compatible with FFT (see Figure [5|) . This is a more detailed understanding of why the meson screening correlators 
are nearly degenerate in Figure [31 

In Figure [6] we show that the correlator C^~^^\z) approaches FFT as the parameter e is tuned to the optimum in 
each of the smearing schemes, approaching closest to FFT at the optimum. We used the distance z = l/T in this 
demonstration because it is neither in the far IR nor in the UV. The optimum HYP and HEX schemes bring the 
correlator closer to FFT than the APE and Stout smearing schemes, although the latter also come very close to FFT. 
Next we explore this difference between schemes. 
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FIG. 6: The correlator C^'^^\z = 1/T) at 2Tc from set N, normalized by its value in FFT, as a function of the smearing 
parameter e in various smearing schemes. 



C. Smearing, taste symmetry and screening masses 
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FIG. 7: The first panel shows screening masses of the local S/PS, V/AV and N at 2Tc as functions of the pion taste splitting 
aSrriTr at 2Tc in data set O. The screening mass of the nucleon has been multiplied by 2/3 in order to compress the vertical 
scale. Each set of screening masses varies linearly with a measure of the pion taste splitting, aSra-,,. The second panel plots 
aSmTT nt T — against the splitting of the corresponding screening masses, aSfips at two different temperatures, but at the 
same lattice spacing; the line y = lOOa;^, is superposed to indicate the slope. 

On examining screening masses, we found that they depend on the smearing parameter e essentially only through 
the taste symmetry breaking measure 

SniTT ~ ni-yr^ji — "^75, (6) 

where the subscripts on the right denote the pion taste structure. The 75 taste is the Goldstone pion. We chose the 
partner with taste structure 757^ as an indicator of taste splitting since it turned out to be relatively easily measured. 
Figure [7] shows the nearly linear dependence of fi/T on SruT^a. The figure shows the clear superiority of the HEX 
scheme over Stout. Using the scaling shown in Figure [7J one could extrapolate screening masses to the limit (5m^ — >■ 0. 
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However, this is premature, since it involves an extrapolation to suboptimal values of e. The S/PS screening masses 
obtained using local operators with dynamical p4 quarks at a comparable temperature turns out to be around 4. ST 

More information can be extracted from the taste-splitting of the screening masses at finite T, 

The only previous study of this kind was reported in @. In Figure [7] we show Sfips as a function of aJm^. In making 
this comparison we held the lattice spacing fixed, with one set of measurements at T = 0, one at T = 2Tc in set O 
and a third at T = 1.33Tc in set P. We find S^ps oc T{a5m-,^Y over the range of values we obtained. This removes 
the ambiguity remarked upon in Q. 

One can argue for this on general grounds. A hadron mass, M, can be written as Ma ~ f^aKjjs^ma^e), where 
we treat e as a generic label for all the parameters which control smearing. A screening mass, can be written as 
^/T = g{aKjj-g,ma,e^ Nt), since Nt = l/{aT), or as afj, = g' (a A-jjs, ma, e, Nt). For data taken at fixed cutoff, aA-^, 
we need not show this parameter explicitly. Although we work at fixed ma, it is profitable to consider the dependence 
on this variable. A series expansion in ma near the chiral limit would yield different dependence for the Goldstone 
pion mass 

am^g = ai\/ma + O (ma)^^'^ , and aAf = 70 + 71 ma + O (ma)^ , (8) 

where M is any other mass scale, and the coefficients depend on e. As a result, aSrriT^ ~ Jo ~ oti^ma + 71 ma and 
5fips/T = 7o + 7jma. If taste symmetry were recovered in the chiral limit by tuning e, then one might argue that 
70 = 7o = and hence aSfips 0: (aSm^)^. This would also mean that all pion tastes are forced to be Goldstones, with 
an expansion starting at order y/ma. Chiral logarithms, which we have neglected here, could become important at 
smaller masses and spoil this scaling. 

Even if smearing achieves a more limited goal of significantly decreasing StUt^ at finite ma without actually recovering 
taste symmetry completely, one might still recover quadratic scaling. All that is needed is that 70 and 7g become 
much smaller than the actual T = taste splitting in the problem. In general one would have 

(5//ps/T - 7o oc (a(5?7i7r - 7o)^- (9) 

The data in Figure [7] shows that 7q and 70, are small compared to aSm^. This quantifies how well smearing works. 
The fact that it seems to work better at finite temperature than at T = with fixed values of aK^rs and ma possibly 
indicates that the Dirac eigenvalue spectrum is simpler. We shall return to this point later. 

The main conclusion from these studies of smearing is the following. Optimising the suppression of UV modes 
automatically improves taste symmetry in the hadron spectrum at T = 0. This leads to superlinear improvement in 
taste symmetry in the hot phase of QCD. In order to gain most from such an improvement, one should then choose 
the best possible smearing scheme. With partial quenching, as here, this would mean working with the optimized 
HYP scheme; with dynamical smeared quarks it would mean working with the optimized HEX scheme. 

IV. RESULTS 
A. Hot QCD: weak coupling and the Dirac spectrum 

We found that the mass-splitting between chiral partners changes rapidly in the low-temperature phase and vanishes 
fairly close to in the hot phase. In Figure[5]we show that A5 and Ay both vanish aXT = LOST^. Also, a comparison 
of sets N and O show very little quark mass dependence at about the same lattice spacing. These results are in contrast 
to the observations in (j, |5| that Ag remains significantly non-zero up to a temperature significantly higher than Tc- 
The change from the old results Q using the same data set P, confirms that the improvement is due to smearing. 

The rapid approach to behaviour similar to weak-coupling theory has implications for the spectrum of the staggered 
Dirac operator. The vanishing of the pion mass in the chiral limit at T = is related to a finite density of the Dirac 
eigenvalues near zero. It was shown earlier in a study of set O with thin-link quarks that a gap developed in the 
massless staggered eigenvalue spectrum a little above Tc, and that the hot phase contained localized Dirac eigenvectors 
p. 

Here we studied the gap by measuring the smallest eigenvalue of the massless staggered Dirac operator, Aq. The 
ensemble average, (Aq), at a given temperature was generally seen to be within a factor of four of the minimum over 
the ensemble. In view of this, we report (Aq). As can be seen in Figure [SI it climbs by two orders of magnitude 
between Tc and 1.06Tc for the smeared Dirac operator. For the thin-link operator, (Aq) rises at significantly higher 
temperature. 
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FIG. 8: Variation of Av^ (left) and As (right) with temperature. Av approaches at the 95% confidence level immediately 
above Tc, independent of smearing scheme and quark mass. 
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FIG. 9: The ensemble averaged smallest eigenvalue of the massless staggered Dirac operator, (Ao), for set N. With optimal 
HYP smearing, the eigenvalue rises by two orders of magnitude in a narrow range above Tc. For the thin- link Dirac operator, 
the rise is much slower. 



One sees some volume dependence in the result. This was studied extensively in [2]| . where it was found that the 
volume dependence becomes negligible when the spatial size, L, is of the order of l/(Ao). For LT = 4, this would be 
at (Ao) — 0.25, which seems to happen at 1.5Tc. In future it would be interesting to study this volume dependence 
further. 

It is also of interest to note that a(Ao) becomes comparable to am at T ~ Tc with optimum HYP smearing. Since 
this happens for all the data sets, it accounts for the lack of quark mass dependence seen in Ay and A5. With 
thin-link quarks this crossing is delayed to T/Tc ~ 1.5, thus affecting all screening phenomena. 

In set N at T = l.STc we spotted one configuration out of the 50 for which Aq was two orders of magnitude below 
(Ao). This implies the existence of a small fraction of atypical configurations in the thermal ensemble. These would 
be interesting in a study of axial U(l) symmetry at finite temperature, where such atypical configurations have been 
linked to topological configurations by observations with overlap [l^l or HISQ quarks in [2^. However, that would 
require a much larger statistical sample, and is therefore best left to the future. 
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FIG. 10: Hadron screening masses for the data sets N, O, and P using optimum HYP smeared correlators. The horizontal lines 
are the free theory screening mass for the nucleon and the mesons respectively: for set N and O they are 5.27 for mesons and 
7.90 for baryons, for set P they are 5.78 for mesons and 8.67 for baryons. DR denotes the weak-coupling prediction of Q|; the 
unlabelled line immediately below this shows the prediction of 0|. 
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B. Comparison with weak-coupling theory 



Finally, the results for the screening masses as functions of T are shown in Figure [10] for all three data sets with 
optimal HYP smearing. Also shown are the values expected in FFT on lattices with the same size. The analysis 
of correlation functions obtained with these smeared valence quarks shows that the screening masses in all channels 
approach FFT at high T. The most striking new feature of this data is that this approach is from above, in conformity 
with the predictions of Q. Similar results are obtained with optimal HEX smeared quarks. We have shown earlier in 
Figure [7] that there is a remaining uncertainty of around 15% in the determination of the the S/PS screening mass. 
This comes from the residual taste symmetry breaking at the best optimization of the screening parameters possible 
at this lattice spacing. Reduction of this uncertainty requires going to finer lattices. 

The weak coupling prediction for the meson-like screening masses is 



Here as is the 2-loop QCD coupling evaluated in the MS scheme at the scale 2itT. Eq — 0.3824 for two flavours 
of quarks in a dimensional reduction (DR) schemed . A hard thermal loop (HTL) resummation which neglects soft 
gluon contributions to the vertex yields Eq = Q . These weak-coupling predictions are also shown in Figure [101 
with as determined using As one can see, both the weak-coupling predictions are close to the observed screening 
masses. 



Several properties of quarks at experimentally accessible temperatures above Tc seem to be explained in weak 
coupling QCD. However, one which showed puzzling departures from weak-coupling predictions was screening masses 
from hadronic excitations. In quenched computations it was seen that the results depended strongly on the kind of 
valence quark used [loj . With this clue in hand we performed computations with dynamical staggered sea quarks and 
improved valence quarks in three sets of computations, one new (set N, see Table [I]) and two older, (sets O [l^ and 
P [13 ) • Studies with staggered valence quarks were reported earlier with set P Q ■ 

A preliminary part of this work was the optimization of the valence quarks. We used four popular versions of 
fat-link staggered quarks. We optimized the smearing parameter, e, in each case by observing changes to the power 
spectrum of the plaquette (see Figure [T|) and the largest and smallest eigenvalues of the Dirac operator (see Figure [2]). 
The optimum e was chosen so that the UV was suppressed as much as possible without changing the IR behaviour 
in both cases. This also improved the performance of the conjugate gradient algorithm used for the inversion of the 
Dirac operator (see Table HI)) . Such a tuning was done at T = 0. We found mild changes in the tuning parameters as 
the lattice spacing was changed by a factor of 2. 

Although the smearing parameter is optimized by requiring that the IR components of fields do not change ap- 
preciably, it docs affect the long-distance properties of the theory, such as masses. We compared different schemes 
through a measure of the recovery of staggered quark taste symmetry in the spectrum of pions (see Figure [T]). The 
optimized HYP smearing works best, although optimized HEX smearing comes a close second. This is pleasant, since 
dynamical simulations with HEX smearing are easier than with HYP. 

Smearing causes systematic changes in finite temperature properties of interest. We found that the screening mass 
in the hot phase increases systematically as taste symmetry breaking is reduced at T = (see Figure [7|). Also, taste 
symmetry breaking in the hot phase improves super-linearly with improvement at T = (see Figure[7]). Since recovery 
of taste symmetry has been used as the main indicator of the reduction of UV effects, it is natural in this study to 
use optimized HYP smearing in order to best reduce lattice artifacts. 

On doing this we find that the screening correlator recovers the degeneracies that a theory of weakly coupled 
fermions would predict (see Figure [3]). This happens very close to, and above, Tc (see Figure [J]). The correlators 
themselves are also close to the predictions of a free fermion field theory (see Figure [5]). Consistent with this, the 
screening masses at high temperature are found to be close to weak-coupling theory (see Figure [TU]). A computation 
in dimensional reduction Q gives results which are slightly different from a HTL computation neglecting soft-gluon 
effects on the vertex Q. The lattice computation is unable to distinguish these as yet, but we may expect this to 
improve in the near future. 

We also see that the smallest eigenvalue of the optimally HYP smeared massless staggered Dirac operator shows 
a rapid jump from extremely small values in the mean below Tc to fairly large values above (see Figure [9]). The 
behaviour of the thin-link staggered operator is qualitatively similar, although quantitatively slower to change. Since 
the smallest eigenvalue of the massless smeared operator is comparable to the bare mass already at T = Tc, the limit 
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FIG. 11: Dependence of rn^xli/T^ on T/Tc- The data points marked by boxes are measurements made on line of constant 
m/Tc, whereas those marked by circles are obtained with constant am. 
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FIG. 12: The spatial size dependence of /3c (left) and A/3c (right). The lines enclose 68% confidence limits on /3c, obtained by 
fitting a single constant to all the estimates. 



of physical renornialized mass becomes easy to take in the high temperature phase. There is evidence for a very small 
fraction of completely atypical configurations in the hot phase. A study of the topology of these gauge configurations 
lies outside the scope of this paper. 

The lattice computations described here were performed on the Cray XI of the ILGTI in TIFR. We thank Ajay 
Salve and Kapil Ghadiali for technical assistance. We also thank Debasish Banerjee, Saumen Datta, Mikko Laine, 
Nilmani Mathur, Subroto Pal and Christian Schmidt for their comments. Some of the configurations used in this 
study were generated earlier for other studies by the ILGTI. 



Appendix A: Determination of /3c 



The cross over is determined at A''^ = 4 for a bare quark mass acjn — 0.015, where is the lattice spacing at jS^- 
We determined Pc by positions of the peaks of different susceptibilities. A/3c was defined to be the full width at half 
maximum (FWHM) of the same susceptibilities. 

We measured the Wilson line susceptibility, Xl the bare chiral susceptibility, Xm [sll, the corresponding 

rcnormalizcd quantity m^x^^/T^ [26l |. and the fourth order QNS, X22 and Xio at various values of /3 in the 
crossover region. For the measurement of m^x^j/T'*, we determined the chiral condensate at zero temperature on 16* 
lattice at the same values of /? as the finite temperature ones. 

Prior to the runs listed in TableUl we performed a series of runs at fixed bare quark mass, am — 0.015, with Ng ~ 8 
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/3c 


A/3c 


Xm 


5.2747(6) 


0.009(2) 


Xl 


5.2743(5) 


0.006(1) 


X22 


5.2741(5) 


0.006(1) 


Xio 


5.2743(6) 


0.007(1) 



TABLE III: Pc and A/3c as determined from different susceptibilities. A/3c is much larger than the statistical error in jic- 



and 12. We used these runs to make first estimates of /3c, and followed up with the runs along lines of constant m/Tc 
listed in Table ID The compatibility of these runs is shown in Figure [TTl where rri^XM/'^'^ is given as a function of 
T/Tc- The figure also shows that with this cutoff, the dcconfining and chiral cross overs in QCD coincide; m^xL/^"^ 
peaks between 0.98 and 1.02Tc. 

To determine Pc accurately, we interpolated data for susceptibilities using multihistogram rcweighting ^3] in the 
cross over region. From bootstrap resampling of the histograms, we determined the means and errors in the position 
of the peak of each susceptibility and its FWHM, so obtaining /3c and A^c [13 • We found /3c and A/3c for each of 
the susceptibilities on the three different lattice volumes. The results of which are shown in Figure [T^l Since we 
found very little volume dependence in /3c, we made a fit to a constant, independent of volume. The values of /3c 
so determined are displayed in Table Hill In Figure [12] we also show the volume dependence of A/3c. This decreases 
with the volume, and gives some indication of saturating, within errors, close to our largest lattice. So wc take A/3c 
obtained on Ng — 16, as our best estimate. These estimates are also listed in Table Hill We find that the variation in /3c 
with different susceptibilities occur well within the width of the cross over measured from each indicator separately. 
In fact, the four estimates of /3c are consistent with each other within 68% confidence limits. Combining all four 
measurements, we quote /3c = 5.2744(7) and A/3c ~ 0.006. 
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